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Abstract 

We set up a formalism to calculate the charge form factors of two-neutron halo nuclei with 

S-wave neutron-core interactions in the framework of the halo effective field theory. The method 

Q-i' is applied to some known and suspected halo nuclei. In particular, we calculate the form factors 

■^ ■ and charge radii relative to the core to leading order in the halo EFT and compare to experiments 

■^ . where they are available. Moreover, we investigate the general dependence of the charge radius on 

the core mass and the one- and two-neutron separation energies. 
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I. INTRODUCTION 

The determination of properties of nuclei along the neutron drip line poses one of the major 
challenges for modern nuclear experiment and theory. The associated observables are an 
important input to studies of stellar evolution and the formation of elements and provide 
insight into fundamental aspects of nuclear structure. Since these systems are weakly-bound, 
drip-line nuclei display novel phenomena associated with newly emerging degrees of freedom 
or, phrased differently, strong correlations. 

Halo nuclei are one example where the transition to new degrees of freedoms becomes 
especially apparent. These nuclei have a tightly bound core with weakly- attached valence 
nucleons \MiM- Usually, they are identified by an extremely large matter radius or a sudden 
decrease in the one- or two-nucleon separation energy along an isotope chain. Thus they 
display a separation of scales which exhibits itself also in low-energy scattering observables 
through a scattering length a that is large compared to the range of the nucleon-nucleus 
interaction R. Halo nuclei can be studied with an effective field theory (EFT) that exploits 
this separation of scales as a small expansionparameter R/a and is formulated in the relevant 
degrees of freedom [a, |6| (See, e.g., Refs. [7|, |8[ for recent reviews.) In this EFT, the core and 
the spectator particles are treated as the fundamental fields in the problem and the overall 
computational complexity decreases significantly. In contrast to ab initio approaches which 
try to predict nuclear observables from a fundamental nucleon-nucleon interaction, halo EFT 
essentially provides relations between different nuclear low-energy observables. On the one 
hand, it thus provides a framework that facilitates a consistent calculation of continuum 
and bound state properties when information on the nucleon-nucleus interaction is known. 
On the other hand, halo EFT can be used to determine two-body scattering properties from 
few-body observables if a sufficient number of them is known. 

In a previous work, Canham and Hammer [9|, [lO(] explored the universal properties and 
structure of 2n halo nuclei to next-to- leading order (NLO) in the expansion in R/\a\ by 
describing the halo nuclei as an effective three-body system consisting of a core and two 
loosely bound valence neutrons. Their main focus was the possibility of such three-body 
systems to display the universal Efimov effect ll| and on the structure of the halo candidates. 



In particular, the matter density form factors and mean square radii were calculated. Using 
this framework, Acharya et al. recently carried out a detailed analysis of the implications of 
a matter-radius measurement for the binding energy and existence of excited Efimov states 
in ^^C [l2|. (For a selection of previous studies of the possibility of the Efimov effect in halo 



nuclei using three-body models, see Refs. |l3Nl6|.) 



Halo EFT has also been extended to include Coulomb effects [17[ and electromagnetic 



currents |18|-|21| in two-body halo nuclei. In this paper, we extend these studies to three- 
body halo nuclei. We consider the electromagnetic charge form factor and charge radius 
of two-neutron halo nuclei interacting through resonant S-wave interactions. While these 
quantities have not been measured yet for current S-wave halo candidates, we anticipate 
that the charge radius will be accessible in the near future due to the constant experimental 
progress in laser spectroscopy techniques on radioactive ions that facilitate a measurement 
of the atomic isotope shift. 

This manuscript is organized out as follows: In Sees. [TTl iHl and HVt we lay out the field 
theoretical formalism required for the calculation of strong interaction observables in two- 
and three-body S-wave halos. The trimer Greens functions are introduced in Sec. |V] and 
the calculation of the charge form factor is discussed in Sec. I VII Our results for the form 



factors and radii including error estimates are presented in Sec. IVIII Finally, we conclude in 
Sec. IVIIII Detailed expressions for the various contributions to the form factor are given in 
the Appendix. 

II. EFFECTIVE LANGRANGIAN 

We set up a non-relativistic effective field theory for a core (c) with spin 0, mass mo and 
electromagnetic charge Ze, interacting with two neutrons (n) with spin 1/2 and mass mi. 
The core is described by a scalar field ipQ and the neutrons are described by a two-component 
spinor field ipi = ( n| ) • 

We assume that all two-particle interactions are short-ranged and dominated by S-wave 
resonances. The treatment of resonant P- and higher partial wave interactions will be left 
for future work. If they are non-resonant, higher-partial wave interactions are suppressed 
by three powers of R/a p, |6|. The core-neutron interaction is described by a spin-1/2 dimer 

field di = i /'^ ] and the interaction of the two neutrons is described by a spin-0 dimer-field 

(Iq. The labeling and indices are chosen to simplify the notation for the three-body equations 
derived below. Moreover, we allow for a three-body contact interaction between the core 
and the two neutrons which is described by a spin-0 trimer auxiliary field t. Note, that 
our choice to introduce auxiliary fields does not imply bound states in the corresponding 
channels and merely is a convenient way to introduce interactions. 

In addition to the strong interactions between the neutrons and the core, we include 
electromagnetic interactions with a vector potential A^. The interaction terms are obtained 
by minimal coupling which insures gauge invariance: idf^ H- idf^ — QA^, where Q is the charge 
operator. In our case only the core has non- vanishing charge Ze, such that QipQ = Zeipo 
and Qipi = holds. For convenience, we choose Coulomb gauge where {V ■ A) = 0. Since we 
restrict our analysis to leading order (LO), non-minimal coupling terms do not contribute. 
The effective Lagrangian can then be written as the sum of one-, two- and three-body 
contributions, C = C^^^ + C^"^^ + C^^\ where 



J, >^ -^ I -^ I >^ , 



(1) 



C^'^ = 4 l^do + ^) ^0 + V^I (*5o + ^ ) ^1 - Ze ^AotjJo 

£(') = Ai didi - gi^dli^iiJo + ijl^fldi 

£(3) =ntU - hlf^ifjodo + iipodoyt] . 

Because we focus on resonant S-wave interactions, the electromagnetic interaction appears 
only in C^^\ 

The two-body part C^'^'' includes the bare dimer propagators and the coupling of a dimer to 
two single particles. The bare parameters Aq, Qq, Ai and gi depend on the ultraviolet cutoff 
A. At LO the parameters Aj and Qi {i = 1,2) are not independent. Physical observables 
only depend on the the combination gf/Ai. The spin projection matrix P projects the 



two neutrons on the spin-singlet. Its components are the corresponding Clebsch- Gordon 
coefficients, leading to 

P = ;^(-io) = -^^ (2) 

such that PP^ = I2/2 and Tr [PPt] = 1. 

Finally, C^^'^ represents the three-body interaction written in terms of a trimer auxilliary 
field J22|. It includes the bare trimer propagator and the coupling of the trimer t to the 
c?o-dimer and the core field ^po■ Writing the three-body interaction using a trimer auxilliary 
field will be convenient for deriving the form factor expressions below. The bare parameters 
fl and h depend on the ultraviolet cutoff A. Again only the combination h'^/Q contributes 
to observables at LO. 

There exists a whole class of equivalent theories in the three particle sector. Integrating 
out the auxiliary fields, one can show that different choices of C^"^^ and C^^^ can be trans- 
formed into the same theory without dimer and trimer fields up to four- and higher-body 
interactions. To demonstrate this, we eliminate the trimer field t using the classical equation 
of motion, resulting in 

£(3) ^ -Hoi^odoYi^odo) , (3) 

where Hq = ^. Repeating this step for the dimer fields di and do then yields 

where Cq = ^-,0^ = ^-, and ifp = ^ ° . The term £'^-^Mncludes interactions of four 
or more particles. In this work, we will only consider processes with at most three particles 
and therefore neglect C'^-^\ By this procedure, physical observables will be unchanged as 
long as the coupling constants are chosen appropriately. In particular, the trimer field could 
also have been introduced in another channel, such as 

£(3) = ^tU - hU^ {^iPdi) + (v?M)^tj , (5) 

without changing any three-body observables. 

In order to write down the derived expressions more compactly, we define the mass 
parameters: 

7. ^ 7. ^ mom? ^ niiMi ,„. 

Mtoi = mo + 2mi , M^ = Mtot - rrii , Hi = — , /ij = — - — . [b) 

rriiMi Mtot 

In the following, we use Feynman rules in momentum space to calculate the properties of 
the en, nn, and cnn systems. In the Feynman diagrams, particles, dimers and trimers 
are denoted by single, double and triple lines, respectively. In addition, propagators are 
represented by arrows, photon couplings by rectangles and all other couplings by ellipses. 
These symbols are empty if they correspond to bare and filled if they correspond to full, 
interacting quantities. Since we consider a non-relativistic theory, the one-body properties 
are not modified by interactions. We thus start with the two-body problem in the next 
section. 
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FIG. 1: Diagrammatic representation of the integral equation for the full dimer propagator. Par- 
ticles and dimers are denoted by single and double lines, respectively. The bare (full) dimer 
propagators are indicated by the empty (filled) arrow. 

III. TWO-BODY PROBLEM 



We start by calculating the full dimer propagator Di{pQ, p) with three- momentum p and 
energy pq for the dimer field with index i = 0,1, corresponding to the nn and en channels, 
respectively. The integral equation is depicted in Fig. [1] in terms of Feynman diagrams. 
Using the Feynman rules derived from the effective Lagrangian ([1]), we find 



iDi{po,p) 
-iT,i{po,p) 



1 z 

— (-iSi(po,p)) iDi{po,p) + ^ = i[Ai-~^i{po,p)] 



■ s^gf^^i 



-1 



■n^ 




(7) 



where s, = 5jo/2 + dn is a symmetry factor and — zSj(j9o,p) is the dimer self energy. It 
is linearly divergent and has to be regularized. For this purpose, we choose a momentum 
cutoff A which is the same in both channels. 

Matching the two-body T-matrix obtained from the propagator ([7]) for p = and po = 
k^/{2fii) to the effective range expansion of the S-wave scattering amplitude 

T 1 -1 



m) 



+ 0{k'^)-ik 



(8) 



with scattering length a,, we eliminate the dependence on the cutoff A. This leads to the 
renormalization condition for the two-particle couplings 

1 27rA,- 2 



-A 



a-i Si gffii 
and the renormalized, full dimer propagator 



TC 



A(po,p) 



2n 



s^ gff^i 



\ 2fii 



P' 



2M 



-Po 



IE 



-1 



(9) 



(10) 



For positive scattering length a,, this propagator has one pole on the first Riemann sheet of 
the complex square root with the positive residue 

Z. = ^,-. (11) 

sigtlii ai 

This pole corresponds to a two-body bound state with binding energy Bi = l/{2fiiaf). For 
negative a, there is a pole with negative residue on the unphysical, second Riemann sheet. 
The leading correction to the propagator flTOj) is due to the effective range. It could be 
included by making the dimer fields dynamical as discussed, e.g., in Refs. 23N25|. Here, we 
stay at leading order in the EFT expansion and neglect effective range corrections. 



IV. THREE-BODY PROBLEM 





+ 




+ 




+ 




FIG. 2: Diagrammatic representation of the integral equation for the particle-dimer T-matrix. 
Particles, dimers, and trimers are denoted by single, double, and triple lines, respectively. The 
bare (full) propagators are indicated by the empty (filled) arrows. 



We proceed by calculating the T-inatrix element Tij for the scattering process of a dimer 
and a particle. The dimer and particle in the incoming (outgoing) channel are both labeled by 
the index i (j), respectively. We focus on the dimer-particle T-matrix as a central quantity, 
since all other three-body observables can be obtained from it. The integral equation for 
the T-matrix is depicted in Fig. [2] in terms of Feynman diagrams. The matrix structure of 
the equation is implicit in Fig. [2j 

The derivation of the corresponding expressions can be performed in several steps. First, 
we project in- and outgoing states onto the spin-singlet using the projection operator from 
Eq. ([2]). External dimer fields are then renormalized with the wave function renormalization 
factors from Eq. fITT]) via 

T,, ^ \Zi\'^Tij\Zj\-^ . (12) 

The absolute values in Eq. flT^ are only required for i,j = because Zq < 0. This channel 
corresponds to a neutron-neutron dimer which is unbound and requires no wave function 
renormalization factor. In this case, Eq. (TT2|) simply provides a convenient redefinition of 
the amplitude but has no physical significance. 

We work in the center-of-mass frame, in which the on-shell T-matrix only depends on 
the total energy E and the relative momenta in the ingoing and outgoing channels p and k, 
respectively. If the dimer and the particle have ingoing (outgoing) three-momenta p^^ and 
P2 (ki and k2) in a general frame, the relative momenta are 



rrii 



M 



M 



-Pi- 



tot 



M 



-P2 



m 



tot 



M 



^ki 



M, 



tot 



M, 



^2 • 



(13) 



tot 



The T-matrix can be decomposed into partial wave contributions 2"[^™'^™]. The Wigner- 
Eckart theorem then implies that T is diagonal in i and m and can be written as T'^'"'^ "^ ^ = 
^w^mm'T^^^. The resulting 2 x 2-matrix integral equation for angular momentum 
generalization of the Skorniakov-Ter-Martyrosian (STM) equation 26| and reads 



IS a 



T^'\E,p,k) = / dqR^'\E,p,q)D{E,q)T^'\E,q,k) + R^'\E,p,k) , (14) 

Jo 



where A is an ultraviolet cutoff on the loop-momentum in the three-particle sector. The 
components of the interaction matrix i?[^l are given through 

pM^p ^ j,\ - (i_^Mz^r ^V ^^ Mtot - m, - m, 1 ^ / .^ ^ ,.^ 
''-^''^P^^^ - ./sis- ^-'^ ^^H ^. -,Q^{c.,{E,p,k)) 

- 5io5jQ 5io H , (15) 

Cij{E,p, k) = I h E — le 



pk \ 2/ij 2/i 

where Qe are the analytically continued Legendre functions of the second kind. In our 
numerical calculations, we will only need 

Qo(c) = /^^^*^^h(l/c) : |c| > 1 , ^^g^ 

1 arctanh(c) + z| : |c| < 1 . 

The form of Qq is determined by taking the limit £ — )■ 0"*^ in the integral equation. Moreover, 
H = \ZQ\h? /Vt is the redefined three-body coupling, which depends on the cutoff A. It only 
contributes for angular momentum £ = 0. The dimer matrix is diagonal in the channel 
indices: D = diag(Z)o, -Di) with 

D,{E,q) = ^— 4===, k{E,q) = 2/i, f ^ - E - zs^ . (17) 



Note that an overall factor — g^/(27r^|Zj|) from the measure of the integration over the loop 
momentum q and the dimer wave function renormalization has been absorbed in Di for 
notational convenience. 

Assuming the existence of an S-wave three-body bound state at energy E = —B, the 
transition amplitude can be decomposed as 

T^'\E,p,k) = -5,0^^^^^-^^ + regular terms . (18) 

E + B + ie 

The residue of the bound state pole factors into wave functions B{p) depending only on 
one single momentum, and the remaining part is a regular function in the energy. Inserting 
Eq. (IT8l) into Eq. ( 1141) yields the bound state equation 

Bip) = [ dq R^'\E,p,q) D{E,q) B{q) . (19) 

Jo 

This generalized eigenvalue problem has an Efimov-like spectrum of three-body bound state 
energies. For a given cutoff A, we then fix the unknown three-body parameter H such that 
Eq. (TT9l) has a solution at the desired value E = —B. In this way, the three-body coupling 
is renormalized and other three-body observables can be predicted. In particular, Eq. ( TT^ 
can be solved numerically in order to determine the T-matrix for three-body scattering 
observables. 

In the following, we will consider only three-body observables in the S-wave {i = 0) 
channel and drop the index "[0]" on the quantities i?'^^ and Tt*^' for notational simplicity. 



From the T-matrix, we can derive the scattering amphtude and scattering length for dimer- 
particle scattering. Since the two-neutron system is not bound, only the element Tn of the 
2x2 T-matrix in Eq. (IT4|) describes a physical scattering process, namely the scattering of 
a neutron from a en bound state: 

\2iii 2/iiaf / III p cot bcn-n[p)-^P 

where the reduced masses /xi and /ii are defined in Eq. ([6]). 

V. TRIMER GREENS FUNCTIONS 
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FIG. 3: Feynman diagrams contributing to the full trimer propagator t. Particles, dimers, and 
trimers are denoted by single, double, and triple lines, respectively. Otherwise the notation is as 
in Fig. [21 

In order to calculate the charge form factor of the three-body bound states corresponding 
to two-neutron halo nuclei, we need the trimer wave function renormalization Z^^ which 
is given by the residue of the bound state pole in the full trimer propagator including 
interactions. The Feynman diagrams contributing to the propagator are shown in Fig. [3l 
For a trimer at rest, it can be written as 

/■A /"A 

t{E,Q) = dq dq'Do{E,q)TooiE,q,q')DoiE,q') + regular terms . (21) 

Jo Jo 

Because the three-body interaction from Eq. ([1]) acts only in the i = channel, only Tqq 
contributes in Eq. (12T!) . The trimer wave function renormalization can be extracted from 
the relation 

Ztr = - lim {E + B) t{E, 0) . (22) 

The trimer self energy is given by all trimer-irreducible contributions to the propagator 
and can be defined as 

E(E) = t{E,0)\^^^, (23) 

where the three-body force is set to zero in the evaluation of Tqq. Using ^{E), the trimer 
propagator can also be written as 

Requiring that t{E, 0) has a bound state pole at £" = —B, Eq. (^^ directly leads to the 
relation 

S(-fi) = 1/H . (25) 

8 
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FIG. 4: Diagrammatic representation of the irreducible trimer-dimer-particle three-point function 
girr_ j'liQ three-body force H is set to zero in the evaluation of the particle-dimer T-matrix. 
Notation is as in Fig. [3l 

In the explicit expressions for the trimer form factor (cf. Appendix |A]), we need the 
irreducible trimer-dimer-particle three-point function defined by 

gr{E,p) := [ dqDo{E,q)Toi{E,q,p) + Soi , (26) 

Jo ^=0 

where again trimer-reducible contributions are removed by setting if = in the integral 
equation for Toj. The corresponding Feynman diagrams are depicted in Fig. |H Inserting 

Eq. ( 1141) into Eq. (1261) and writing Q^^^ = ( "„ j , yields a matrix integral equation for the 

irreducible trimer-dimer-particle three-point function 

g^'^'^iE^p) = [ dqR{E,p,q) D{E,q)g'"{E,q) + e^ , (27) 

Jo ^=^ 

where (eo)j = 5oi. ^'" depends on the cutoff A, but the combination -v/Zt^if ^^" is indepen- 
dent of A up to an overall sign. This is exactly the combination that enters into the form 
factor calculation. 



VI. CHARGE FORMFACTORS 

We are now in a position to calculate the charge formfactor J-e of a cnn halo system with 
resonant S-wave interactions. The form factor can be extracted from the matrix element of 
the electromagnetic current j^ between trimer states. We will denote the in- and outgoing 
three-momenta of the trimer by P and K, respectively. It is convenient to extract the charge 
form factor from the matrix element of the zeroth component of the electromagnetic current 

(t(iro,K)|jo|t(Po,P)) = (-^eZ) J-e(Q') (28) 

in the Breit frame, where no energy is transferred by the photon, i.e. Pq = Kq and P = K . 



The charge formfactor in Eq. (l28l) . then depends only on the three-momentum transfer 
Q2 = (K-P)2. 

The LSZ reduction formula implies that the current matrix element in Eq. (|28|) can be 
expressed as J27l | 

(29) 



:t(Ko,K)|jo|t(Po,P)) = VZtrS(-fi)-l ^ro(Q) S(-fi)-V^tr 



where iro(Q) is the sum of all irreducible Feynman diagrams with external trimer lines and 
a photon coupled to the core. 



FIG. 5: Exemplary irreducible graph contributing to the trimer form factor. The dashed line 
represents the core field which either (a) propagates parallel to the do-dimer, (b) appears within a 
di-dimer loop or (c) is exchanged between di-dimers. Otherwise the notation is as in Fig. HI 

In order to motivate the different contributions to zro(Q), we consider the typical irre- 
ducible graph shown in Fig. |5l The photon can only couple to the core field c indicated 
by dashed lines, but for the moment we suppress the photon-core coupling.^ Within such 
a diagram, the core (a) propagates either parallel to the do-dimer, (b) appears within a 
c?i-dimer loop, or (c) is exchanged between between two (ii-dimers. In fact, these are the 
only 3 possibilities for a single core propagator to appear in an arbitrary irreducible trimer 
graph. Thus, including the photon-core coupling and summing over all such diagrams the 
form factor derived from the transition amplitude through Eq. f l28l) and Eq. ( I29l) . can be 
written as the sum of three contributions 



^F 



^1 



(a) 






^ -/ E 



(30) 



corresponding to the cases (a), (b), and (c). In all three contributions the irreducible trimer- 
dimer-particle three-point function Q^" from Eq. (126|1 appears naturally. 





L 



(a) 





FIG. 6: Diagrammatical representation of the form factor matrix element (j28p . The contributions 
fall into three different classes (a), (b), and (c). Notation is as in Fig. SI 



In Fig. El the decomposition of the form factor matrix element ( l28ll into the three classes of 
diagrams is illustrated pictorially. The filled circle represents Q"'^ from Eq. ( !26l) . Performing 
shifts in the loop momenta of order Q = |Q| <C A one can express J^ ' *^e ' ^^^ -^e 



^ Of course the di-dimer also carries charge, but the photon coupling to di appears only at next-to-leading 
order where the dimers are dynamical. 
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through integrals symmetric in the momenta. The various contributions are derived in 
Appendix [X] and exphcit expressions are given in Eqs. ( 1A3|) . (1A6|) and (1A8|) . respectively. 

At Q^ = 0, the charge formf actor is normalized to one. This normalization is automat- 
ically reproduced in our formalism. For small momentum transfer, the form factor can be 
expanded in powers of Q^ as: 

MQ') = 1 - ^Q' + o{Q') , (31) 

where (r^) is the charge radius. 

In practice, we calculate the form factors for finite momentum transfer and extract the 
the charge radius (r|) by numerically taking the limit 

(rl) = -6 hm ^ . (32) 

Here, we have to keep in mind, that in our effective theory the core and the neutrons are 
treated as pointlike. Their size enters only in counter terms that appear at higher orders. 
In typical halo nuclei, however, the charge radius of the core can not be neglected. In this 
work, we thus interpret the calculated radius as the charge radius of the cnn halo nucleus 
relative to the charge radius of the core. The small negative charge radius of the neutron 



(r|)^ = —0.115(4) fm^ 28| is neglected. In order to get the full charge radius of the cnn 
halo nucleus (r|)^^^, we therefore have to quadratically add our result to the charge radius 
of the core (r|) : 

(-•D^ = (r|>, + (r|> , (33) 

or simply quote the difference 5 {r"^) := {r^)cnn ~ (^e)c- '^^^^ prescription follows directly 
if the total charge distribution is a convolution of the charge distributions of the halo and 
the core. The differences of nuclear charge radii 6 (r|) were measured for a whole range of 
isotopes, see e.g. |29|-|33l|. In the next section, we will compare these results with our theory 
where it is applicable. 



VII. RESULTS 

We now apply our effective field theory to concrete physical systems. Our theory applies 
directly to two-neutron halo nuclei with J'^ = 0"^ and with a J'^ = 0~^ core. Assuming that 
the spin of the core is inert due to the large mass of the core compared to the neutrons, we 
can also consider more general systems with quantum numbers J'^ , (J ± 1/2)^^, and J'^ of 
the C-, en-, and cnra-systems, respectively. From now on, we write m„ for the neutron mass 
and rric for the core mass for convenience. 

In order to have a bound or virtual cn-dimer and a bound cnn-trimer their one- and 
two-neutron separation energies Ben = Bi and Bcnn = B have to obey Bcnn > niax(Scn, 0). 
Moreover, we denote the first excitation energy of the core by E* and the one-neutron 
separation energy by Bc-n- As discussed in the introduction, the expansion parameter of 
our theory Rcore/Rhaio is roughly R/a. In order to obtain better estimates, we compare 
the typical energy scales -Ehaio and -Ecore of the neutron halo and the core, respectively. 
To estimate -Ehaio, we choose the one- or two-neutron separation energy B^n or Bcnn- The 
energy scale of the core is estimated by its excitation energy E* or its one-neutron separation 

11 



energy -Bc-n- The square root of the energy ratio -Rcore/-Rhaio ~ A/|-E'haio/-E'core| then yields 
an estimate for the expansion parameter of the effective theory. For our error estimates, we 
take the largest value for -RcoreZ-Rhaio that can be obtained this way. 

We fix the values of all quantum numbers J\ masses, and energies by taking data from 
the National Nuclear Data Center (NNDC) [35|. The cn-scattering length is determined 
from the relation Oi = a^n = sign(i?i)/A/2/ii|i?i|. Thus cn-states with negative Bi are 
treated as virtual two-body states with negative scattering length ai. This approximation 
corresponds to neglecting the imaginary part of the binding momentum for resonances. For 
the ?7,n-scattering length, we take the value a^ = ann = —18.7(6) fm from Gonzales Trotter 
et al. Q. 

The lightest isotopes for which there is either experimental evidence for their 2?T,-halo 
nature or which are good candidates for such a system, are ®He, ^^Li, ^^Be, ^^B and ^^C 
Since the J'^-quantum numbers of ^He and ^^B indicate that P-wave contributions must be 
dominant, we apply our effective theory to ^^Li, ^^Be and ^^C 

In Tab. [U we summarize the effective theory parameters and our predictions for the 
charge radii relative to the core in ^^Li, ^^Be and ^^C Uncertainties in the energies are only 
quoted if they are larger than 1%. The small uncertainties in the nuclear masses can be 
neglected. For ^^Li and ^'^Be, the expansion parameter Rcore/Rhaio is typically not much 
smaller than 1. As a consequence, the main uncertainty in our calculation for these systems 
is from the next-to-leading order corrections in the effective theory. In the ^^C halo nucleus, 
the main uncertainty is from the poor knowledge of the binding energies. Below, we discuss 
our analysis for each halo nucleus in detail. 
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3" 
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0.11(6) 





2.9(3) 

-0.072 
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l-66to~49 



TABLE I: Effective theory parameters, estimates of the expansion parameter, and predicted electric 
charge radii relative to the core 5 {r^) from Eq. ([32|) for the halo nuclei ^^Li, ^^Be and 22c. Further 
explanations are given in the text. 



The 



while 



^^Li halo nucleus and the ^Li core have both the quantum numbers J'^ = | 
^°Li appears to have either J'" = 2~ or 1~. There is some evidence that both S- and P- 
wave components contribute to the neutron halo [2|. However, we will analyze ^^Li under 
the assumption that only the S-wave contributes in LO and test the consistency of our 
assumption with the data. P-wave contributions will enter in higher orders. The two- 
neutron separation energy of ^^Li is 0.37 MeV and ^''Li is 26(13) keV above the n-^Li 
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FIG. 7: (Color online) The charge form factor J^e(Q^) for the halo nuclei ^^Li (left plot, blue solid 
line), ^^Be (left plot, red long-dashed line) and '^'^C (right plot, green dash dotted line) relative to 
the core in leading order halo EFT. The estimated theory error for ^^Li and ^^Be is given by the 
shaded bands. For '^'^C, varying the separation energies within their errors gives the shaded region. 
The vertical dashed lines indicate the breakdown scale from explicit pion exchange. 



threshold. The first excitation energy of the ^Li ground state is 2.69 MeV and its one- 
neutron separation energy is 4.06 MeV. Thus the expansion parameter and the error can 
be estimated as Rcore/Rhaio ~ \/Bcnn/El ~ 0.37. Calculating the charge radius relative to 
^Li via Eq. fl32|) gives 5 (r|) = 1.68(62) fm^, where the ~ 40% uncertainty comes from the 
expansion parameter. In Ref. [291, the charge radius was measured with the help of high 
precision laser spectroscopy. The experimental value of 5 {f'^ = 1.171(120) fm^ is thus 
compatible with our calculation within the error bars. 

The halo nucleus ^^Be and its core ^^Be are both in a J'^ = 0+ configuration, while the 
quantum numbers of ^^Be are less clear although there is some evidence for J'^ = ^ . For 
our study, we assume that the ^^Be dimer has also positive parity. The binding energy of 
the ^^Be trimer is B^nn = 1.27(13) MeV and the virtual ^^Be has B^^ = -510(10) keV. The 
excitation energy of the ^°Be core is E* = 2.10 MeV and its one-neutron separation energy 
is 3.17 MeV. Thus, the resulting expansion parameter Rcore/ Rh&\o ~ \/ BcnrJ E^ ^ 0.78 is 
relatively large. Using Eq. f l32|) . our effective theory then predicts a charge radius relative 
to i^Be of 5 (r|) = 0.41(32) fm^ with an ~ 80% error. 

There is some theoretical and experimental evidence that ^"^C is a pure S-wave halo 

= 0+, while 2iC is in J^ = |+ 

0.11(6) MeV has a relatively large 

= -0.014(467) keV is only poorly 



nucleus 36, 37 



2^0 and the ^^C core both have J" 
configuration. The two-neutron separation energy Bcnn - 
error. Furthermore, ^^C seems to be unbound, but 5™ 



13 



known. Up to now, no excitations of the core were observed. The one-neutron separation 
energy of ^°C is 2.9(3) MeV. We take the central values for Bcnn and Ben, which are also 
roughly in accord with the allowed parameter region predicted from a recent analysis of the 



37 



matter radius measurement 
charge radius relative to ^°C via Eq. 



38(1 . Calculating the 



in the framework of the halo EFT 

!]) gives 6 (r|) = 1.661^49 ^^y where the uncertainty 
now comes from varying the separation energies within their errors. Due to the poorly known 
input data, Ben = B^nn is not excluded. Since for such values, the charge radius diverges 
towards positive infinity, the predicted value for 6 (r|) can only be bound from below, where 
the lower limit is 1.17 fm^. For the halo nuclei ^^Be and ^^C, our results are true predictions 
and can be compared with measurements as soon as the corresponding experimental data 
is available. 

In Fig. [71 the charge form factors calculated from Eq. fl30|) are depicted as a function of 



the momentum transfer Q . For Q -^ oo, the form factors vanish. At small momentum 
transfers, they approach unity as required by current conservation. Numerical deviations 
from unity at vanishing Q^ are less than 10~^. This provides a consistency-check for our 
calculation. For ^^Li and ^^Be, the estimated error from higher orders in the effective theory 
expansion is given by the shaded bands. For ^^C the shaded region originates from varying 
the binding energies within their errors. Our effective theory does neither include explicit 
pion dynamics nor does it include the structure of the core. Thus it breaks breaks down 
for momentum transfers of the order of the pion mass m^ ^ 0.5 fm~^ as indicated by the 
vertical dashed lines in Fig. [71 
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FIG. 8: (Color online) The electric charge radius relative to the core 5 (r'^ as a function of mass 
ratio rric/ran for different binding energies Ben and Bcnn- 



Another interesting aspect is to understand the general dependencies of the charge radius 
on the core mass and the en and cnn binding energies. In Fig. [HI we therefore plot 5 (r|) 
as a function of the mass ratio rric/mn for fixed energy values B^nn and B^n- The charge 
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FIG. 9: (Color online) The electric charge radius relative to the core (5 (r^) as a function of the 
energies Bcnn and 1 — Bcn/Bcnn for fixed mass ratio mc/mn = 10. 

radius is always positive and decreases for growing core mass rric. This reflects that 6 (r|) 
for the two- neutron halo system originates from the recoil effect of the charged core. For 
core masses below 2m„ the radius roughly falls of like ~ l/rric- Around rric ^ 2m„, the slope 
changes and the charge radii drop even faster as m,c/mn increases. 

In Fig. [HI the dependence of the charge radius relative to the core on the energies Bcnn and 
1 — Bcn/Bcnn IS showu for a fixed mass ratio rric/mn = 10. As one would also naively expect, 
6 (r|) grows as both the binding energies for the three-body system Bcnn and the binding 
energy of the dimer-particle system Bcnn — Ben decrease. However, the exact functional 
dependencies on the 3 quantities tjIc, Bcnn and Ben are more complicated. Also note that 
there is a sudden increase in 6 (r|) along the line Ben = where the cn-system changes from 
a bound state to a virtual state. This leads to a ridge along {Bcnn — Ben) /Bcnn = 1 that is 
most easily seen in the contour lines. 



VIII. SUMMARY AND CONCLUSION 

In this paper, we have set up a formalism to describe the electromagnetic structure of 2n 
halo nuclei within halo EFT. We assumed that the en and nn S-wave scattering lengths are 
much larger than the range of the interaction and calculated the form factors and charge radii 
of various halo nuclei to leading order in the expansion in Rcore/Rhaio- The renormalization 
of the trimer propagator and the extraction of the trimer wave function renormalization were 
discussed in detail. The charge form factor receives contributions from three different classes 
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of diagrams illustrated in Fig. [61 In all three contributions, the irreducible trimer-dimer- 
particle three-point function Q""^ from Eq. (126!) appears naturally. Current conservation 
insures the correct normalization of the charge form factor J-e(O) = 1. Numerically, we find 
deviations from unity at vanishing Q^ of less than 10~^ in our calculation. 

We have applied our formalism to ^^Li, ^^Be and ^^C calculated their charge form factors 
and radii relative to the core to leading order in -RcoreZ-Rhaio- The resulting charge radii are 
S (r|)nLi = 1.68(62) fm^, S (r|)i4Be = 0.41(32) fm^, and 6 {rl)22c = l-66;;^49 f^^. For "Li a 
comparison with the measured value 1.171(120) fm^ shows good agreement within the ~ 40% 
uncertainty originating from the expansion parameter of our leading order calculation. The 
other charge radii are true predictions that can be compared to future experiments. For 
a more quantitative comparison with experiment, the extension to higher orders is clearly 



required. This includes the treatment of effective range effects |23l-l25| as well as of P- 
wave interactions in ^^Li and ^^Be. Electromagnetic breakup reactions can reveal additional 
information on the structure of 2n halo nuclei. An investigation of this process in the 
framework of halo EFT is in progress [39J] . 

Finally, we have investigated the dependence of the charge radius on the core mass and 
the cnn and en binding energies. Our results are summarized in Figs. HI and |9l As expected, 
the charge radius decreases with increasing core radius. However, the exact dependence for 
for large core masses deviates from a simple l/rric dependence. Moreover, the charge radius 
increases as binding energies Bcnn — Bcn and Bcnn decrease. In particular, we found a sudden 
increase of the charge radius along the line Ben = where the cn-system changes from a 
bound state to a virtual state. A better understanding of these characteristics will require 
further studies. 
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Appendix A: Form factor contributions 

In this appendix, we explicitly derive the three different contributions to the charge form 
factor discussed in Sec. |Vll Therefore, we first give an expression for 2ro(Q) from Eq. (12^ 
as the sum of all irreducible Feynman diagrams with external trimer lines and a photon 
coupled to the core. It can be conveniently written as 

^ro(Q) 



(2vr)^ J (27r)4 .^ . 

lpl<A |k|<A y^^i 

X z^"*-(E,P,po,p)^^fo(S,P,Po,P,K,fco,k) z^'"(E,K,fco,k) . 

The quantity Qf^{E, P,po; p) is the irreducible trimer-dimer-particle three-point function in 
general kinematics, where the trimer, dimer and particle four-momenta are written as P^, 
■^P'' +P^ and -^^P^ - P^, respectively. The energy reads E = Pq - P^/(2Mtot), where 
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the kinetic energy of the three-body bound state is subtracted. Qf^{E, P,po, p) is related to 
the center-of-mass quantity Qf^^{E,p) from Eq. (126|) through the integral equation 



L /-A . /l/T -D . ^ ^2 



r / M P n "n 



H=0 



(A2) 



X D,iE,q)gfiE,q) + 5, 



OJ 



The matrix valued function iTq in Eq. flAll) is the sum of the three diagrams that are 
depicted in the lower row in Fig. [61 Its contributions (a), (b), and (c) are products of the 
corresponding delta-functions, propagators and couplings. 

Since in the Breit frame in- and outgoing three-body bound states have E = —B, we 
drop this redundant energy variable in ^'" and D in all subsequent equations. 

1. Contribution T^"'^ 



(a) 





FIG. 10: Type-(a) contribution to the form factor matrix element where the photon couples to a 
core field propagating parallel to the do-dimer. Notation is as in Fig. HI 

We start with the type- (a) contribution depicted in Fig. [TO] where the photon couples 
to a core field propagating parallel to a c?o-dimer. All non-perturbative physics is con- 
tained in the irreducible trimer-dimer-particle three-point function Q^". One of the two 
four-momentum integrations in Eq. (lAip is absorbed by a delta-function and the remaining 
loop four-momentum q^ can be chosen in such a way that the resulting expression is sym- 
metric under q i— ;■ — q. The g° integration can then be performed analytically, leading to two 
contributions from picking the poles of the two core propagators. The trimer-dimer-particle 
three-point functions Q"'^ appear in off-shell kinematics and are related to on-shell center- 
of-mass quantities via Eq. flA2p at the expense of two additional momentum integrations. 
Choosing spherical coordinates for q, with Q = K — P pointing in the z-direction, the az- 
imuthal dependence can also be integrated out. We are then left with a remaining integral 
over q = |q| and the polar angle x = cos(Z(Q, q)). Since the do dimer is propagating parallel 
to the core, the integral equation f[27|) leads to extra contributions from the three-body force. 
Finally, we end up with 



+ 2J\pg^"{pfD{p)T^'^{Q,p) + Tt\Q) 



(A3) 
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where the matrix-, vector- and scalar- valued functions T^"'\Q,p, k), T^"'\Q,p) and Tq (Q), 
using the shortened notation T*^")(Q . . . ), are given by: 



t^^^Xq. 



4 Jo q 7-1 X ^- VMtot2'^'''' 



(A4) 
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(A5) 



with d{r,q,x) = a/q'^ + 2qrx + r^. One can easily check that the symmetry T^"'\Q,p,k) = 
T^°-\Q,k,p)'^ holds as required. Note also that the limit Q — > in Eq. (JASP exists, but 
prefactors oc 1/Q cause numerical instabilities for very small momentum transfer. 



2. Contribution J"W 



(b) 




FIG. 11: Type-(b) contribution to the form factor matrix element where the photon couples to a 
core field inside a nc bubble. Notation is as in Fig. [H 

We now consider contributions to the formfactor of type (b), where the photon couples to 
a core field inside a nc bubble. The corresponding Feynman diagram is depicted in Fig. [TTJ 
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The bubble sub-diagram can be calculated analytically using Feynman integrals. Analogue 
to case (a) one of the two four-momentum integrations in Eq. (jAip is absorbed by a delta- 
function and the remaining loop four-momentum q^ can be chosen in a symmetric way. 
Applying Eq. (]A2p then again leads to two additional momentum integrations. In spherical 
coordinates, the azimuthal integration can then be performed leading to two remaining 
integrals over q = |q| and the polar angle x = cos(Z(Q, q)). Since diagrams of type (b) only 
appear in channel 1 where the neutron-core dimer is present, there are no three-body force 
contributions in the integral. We find 



fA rA 



J^^'\Q^) = Z,,H^ [ dp [ dkg'"{pf D{p)T('\Q,p,k)D{k)g'''{k) , (A6) 

Jo Jo 

where the matrix-valued function T^^\Q,p, k) is given through: 



Rtii-B,p,d{s,q,-x)) 






arctan ' ^ ^ 



MuA mo_ \^ Ri.j(—B,d(s,q,x),k) 

-|- arctan 



H=0 



y^bi{-B,d{s,q,x))J j -^ + \/bii-B,d{s,q,x)) 

Again T^^\Q,p, k) = T^''\Q, k,p)^ holds and the limit Q ^ can lead to numerical insta- 
bilities. 



3. Contribution J"(^) 




FIG. 12: Type-(c) contribution to the form factor matrix element where the photon couples to a 
core field exchanged between di dimers. Notation is as in Fig. HI 

In the remaining contribution of type (c) shown in Fig. [121 the photon couples to a 
core field that is exchanged between di dimers. For both loops, the energy integrals can 
be performed analytically, leading to on-shell conditions for the trimer-dimer-particle three- 
point function ^"'''. In spherical coordinates, one of the two azimuthal integrals can be solved 
analytically, such that in the end five integrals remain. For J-'^^\ we then get 

^(^)(g2) = z,,H^ f dp f dkg^'^ip)'^ D{p)T'^'\Q,p,k)D{k)g"'{k) , (as) 

Jo Jo 
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with a the matrix- valued function T^^\Q,p, k) given by: 



Tg)(Q,p,fc) =^ f dx f dy r 
\ai\ J_i J_i Jo 



# x\f f -j^Q' ^' y^ <^' P' ^ 



x[f{t,x,y,(f),p,k) = 6ii6ij 
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+ 2ixiB 
+ 2fiiB 



(A9) 



As in the other cases, the symmetry condition T'^^^(Q,p, fc) = T'^^^(Q, k,p)^ holds. 

The dependence of all three formfactor contributions (a), (b), and (c) on Q^ only is not 
directly evident from Eqs. (IA3p . (IA6p and (lASp but numerically it is satisfied. Moreover, the 
charge form factor J^(Q^) is automatically normalized to unity at zero momentum transfer. 
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